ITP-SB-93-63 
November, 1993 



General Determination of Phases in 
Leptonic Mass Matrices 



Alexander Kusenkcj] and Robert ShrockQ 

Institute for Theoretical Physics 
State University of New York 
Stony Brook, N. Y. 11794-3840 



Abstract 

We construct new invariants and give several theorems which determine in general (i) the 
number of physically meaningful phases in charged lepton and neutrino mass matrices and (ii) 
which elements of these matrices can be rendered real by rephasings. We illustrate our results with 
simple models. 
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Understanding fermion masses and mixing remains one of the most important outstanding prob- 
lems in particle physics. In particular, the issue of possible neutrino masses and associated lepton 
mixing is of fundamental interest. Although there is no definite direct evidence for nonzero neutrino 
masses)^ they are expected on general grounds: given only the known left-handed neutrino fields 
and the usual Higgs field(s)PL nonzero neutrino masses result generically from higher-dimension 
operators which would occur at a scale near to that of quantum gravity, suppressed by associated 
inverse powers of the (reduced) Planck mass, Mp = \fhcj (8ttGn) = 2.4 x 10 18 GeV, For exam- 
ple, gauge-invariant dimension-5 operators of this type could produce neutrino masses of order 
m u ~ av 2 /Mp, where v = 250 GeV is the scale of electroweak symmetry breaking, and a is a 
dimensionless constant. Thus, one can understand on general grounds why neutrino masses are 
so small. If electroweak-singlet (EWS) neutrino fields exist, they would lead, via renormalizable, 
dimension-4 operators, to neutrino masses m u ~ v 2 /Mp, where the scale Mp of the EWS neutrino 
mass is naturally » v, again yielding, albeit for a different reason, very small m u [0]. In turn, 
a natural concomitant of (nondegenerate) neutrino masses is lepton mixing, which is thus also a 
generic expectation.^ It is thus of interest to study models for the leptonic mass matrices in the 
charge Q = and Q = — 1 sectors. The diagonalization of these mass matrices determines both 
the masses and the observable lepton flavor mixing. In analyzing models of fermion masses and 
mixing, an important step is to determine the number of real amplitudes and unremovable, and 
hence physically meaningful, phases in the mass matrices. Surprisingly, there was no general solu- 
tion in the literature to the problem of counting the number of complex phases and determining 
which elements of mass matrices can be rendered real by fermion rephasings. We recently presented 
a general solution for the quark sector The leptonic case is qualitatively different and more 
complicated because of the general presence of both Majorana and Dirac masses. In this Letter we 
give a solution for the leptonic case. 

The mass terms for the charged leptons can be written in terms of the Gsm = SU(3) x SU(2) 
x U(l) fields as 

3 

\L 2j LM^£ kR \ + h.c. (1) 



I"/ i r il] , 

j,k=l 



where j, k = 1,2,3 denote generation indices, Ljp = ( j ) , l\ = e, 1% = [i, £3 = r, the first 



1 For reviews and current limits, see Refs. [jij and The apparent solar neutrino deficit is the most suggestive 
indirect evidence at present. The situation with atmospheric neutrinos is unclear. 

2 This comment applies formally to the standard model; however, the Higgs sector of this model is presumably 
stabilized by supersymmetry, in which there are (at least) a pair of Higgs fields with Y = 1 and Y = — 1. 

3 The lepton mixing angles are functions of ratios of elements of neutrino matrix elements and of charged lepton 
mass matrix elements, and even though left-handed neutrino masses are small, some of these ratios could, in principle, 
be O(l). However, a set of conditions for natural suppression of observable lepton flavor violation were formulated, 
and it was shown that the standard model (generalized to include nonzero m v ) satisfies these [ffl. 
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+ h.c. (2) 



subscript on L a jL is the SU(2) index and is the charged lepton mass matrix.^] In contrast 
to the charged lepton sector, where one at least knows the relevant fields, in the neutral lepton 
sector, one does not. In addition to the three known left-handed / = 1/2, I3 = 1/2 Weyl neutrino 
fields VjL, j = 1 — 3, there could be some number n s of electroweak-singlet neutrino fields Xj,R, 
j = 1, 2, ...n s . The general neutrino mass matrix is given by 

(M(°)) T 

where v h = (y e , VtYl, X r = ixii ■■■■>Xn B ) 1 R, Ml and Mr are 3x3 and n s x n s Majorana mass 
matrices, and Mp is an 3-row by n s -column Dirac mass matrix^] We denote the (3 + n s ) x (3 + n s ) 
neutrino mass matrix in (|) as M<°). Thus, M$ = m]° } , j, k = 1,2,3; Afffi = Af£J i>3+fcJ 
j,k = 1, ...n s ; and = M^- k , j = 1, 2, 3, k = l,..,n s , . In general, all of these matrices are 

complex. Recall that the anticommutativity of fermion fields in the path integral and the property 
C T = — C (where C is the Dirac charge conjugation matrix) imply that 

M U) = M (f)T ^ j Qr f = l,R (3) 

so that 

M (0) = M (0)T (4) 

The diagonalization of (§) yields, in general, 3 + n s nondegenerate Majorana neutrino mass eigen- 
states (some could, of course, be degenerate in magnitude, leading to possible Dirac neutrino 
masses). 

Although our theorems on phases are general, some comments are necessary concerning possible 
electroweak-singlet neutrinos. In the supersymmetric generalization of the standard model, which is 
strongly motivated by its ability to maintain the condition v « Mp beyond tree level without fine- 
tuning, the component fields x'j l wou ld arise as part of G^M-singlet chiral superfields Xj ( an chiral 
superfields will be written as left-handed). One may recall that, in general, chiral superfields which 
are singlets under the standard model gauge group Gsm can destabilize the hierarchy Q . However, 
in commonly used models, the symmetry (such as matter parity) which prevents excessively rapid 
proton decay also excludes the types of terms which would destabilize the hierarchy. 

To count the number of physically meaningful phases, we rephase the lepton fields so as to 
remove all possible phases in and by 

L jL = e-^L' (5) 



4 Here we use the result from LEP and SLC that there are three generations of fermions with associated light 
neutrinos Jf). 

5 In the context of the (not necessarily minimal) supersymmetric standard model, we assume unbroken R parity 
so that the neutrinos do not mix with the neutralinos (higgsinos and neutral color-singlet gauginos). 
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t jR = e^i'n (6) 



for j = 1, 2, 3, and, if there exist any XjR s > a l so 



XjR = e J Xjr (7) 



for j = 1, ...n s . In terms of the primed (rephased) lepton fields, the mass matrices have elements 

M f = e^ + ^M$ (8) 
for the charged leptons, and, for the neutrino sector, 



M (L), = e i(a J+Qk ) M (L) (g) 



for j,k = 1,2,3; 



for j, k = 1, n s , and 



M (R)> = ji^+^Mf (10) 



for j = 1, 2, 3 and fc = 1, n s . 

Thus, if AfW has N e nonzero (and, in general, complex) elements, then the equations for 
making these elements real are 

aj + $ = -arg{M$) + r^vr (12) 

where the set {jk} runs over each of these nonzero elements, and ry- k = or lj^| 

Similarly, in the neutrino sector, if M^ D \ M^ L \ and M^ R ' have, respectively, Nd, Nl, and Nr 
nonzero (and, in general, complex) elements, then the corresponding equations for making these 
elements real are 

a j + /3% ) =-arg(M£ ) )+ V $K (13) 
a . +ak = -arg(M^) + t/^tt (14) 

and 

+ /#> = -arg(M^) + r,$\ (15) 

where the various ranges of indices are obvious from (|2|). Let us define the (6 + n s )-dimensional 
vector of fermion field phases 

v = ({<*}, {ti%{ti x) }) T (16) 

6 The rjjk term allows for the possibility of making the rephased element real and negative; this will not affect the 
counting of unremovable phases. 
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where {on} = { ai ,a 2 ,a 3 }, {/f >} = and {/^ x) } = •«,/$?}■ We also define 

the vector of phases of elements of the various mass matrices 

w = ({arg(M^) + t^tt}, {{an^ )) + r?j° fe ) vr}) T (17) 

of dimension equal to the number of rephasing equations N eq = N^+Nq, where Nq = Nl+Nr+Nd. 
We can then write fll2|)- (15) as 

Tv = w (18) 

which defines the N eq -iow by (6 + n^-column matrix T. 

Our first main theorem is: The number of unremovable and hence physically meaningful phases 
N p in AfW and M<°) is 

iVp = N eq - rank(T) (19) 

In order for these phases to be physically meaningful, it is, of course, necessary that they be 
unchanged under the full set of rephasings of fermion fields. In fact, as we will show, there is a one- 
to-one correspondence between each such unremovable phase and an independent phase of a certain 
product of elements of mass matrices which is invariant under fermion field rephasings. (In special 
cases a model may have an unremovable, invariant phase which, because of a particular symmetry, 
is zero or tt; we give an example below.) The proof of ( |l9|) is similar to that for the analogous 
theorem which we proved for the quark sector ||: Let rank(T) = r T . Then one can delete N eq — r T 
rows from the matrix T, i.e. not attempt to remove the phases from the corresponding elements 
of the AfW and M<°). For the remaining r T equations, one moves a subset of 6 + n s — (N eq — r T ) 



phases in v to the right-hand sides of (jl2|)-(15), thus including them in a redefined w. This yields 
a set of r T linear equations in r T unknown phases, denoted v. We write this as Tv = id. Since, by 
construction, rank(T) = r T , T is invertible, so that one can now solve for the r T fermion rephasings 
in v which render r T of the N eq complex elements real. Hence there are N eq — r T remaining phases 
in AfW and as claimed. □{] 

Some comments are in order. First, as is clear from our proof, in general, the result ( |l9D does 
not depend on whether or not = Mff initially for f = £,D. Hence making or 

initially (complex) symmetric does not, in general, result in any reduction in N p , since in terms of 

the Gsm theory, nothing guarantees this symmetry, and the rephasing is carried out in terms of 

. (f) 
Gsm fields. Second, if one of the unremovable phases is put in a given off-diagonal M^ q > for / = t 



7 In the quark case, T is a N eq x 9 matrix, and we noted |]] that it could not have maximal rank because one 
overall rephasing left the Yukawa (or equivalently, mass) matrices invariant. In the present case, unless and 
M^ RS) both vanish, there is no analogous overall rephasing which leaves and M 1 - ' invariant, and hence T can 

have maximal rank. 
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or D, one may wish to modify the qp'th equation to read 

<Xq + 4 f) = -arg(M$ ) - arg(M$ ) (20) 
for f = £,D. For example, in a model where \MpP \ = \MqP\, this would yield Mpp* = Mqp for this 



pair pq and / = £ or D. The modification in (|20|) has no effect on the counting of phases. Of course 
rty 



given the property (|3|), when one rephases to make Mjl real for f = L,R, this automatically does 



the same for M k - . 

A fundamental question concerns which elements of and (taken to be complex ini- 

tially) can be made real by fermion rephasings. This is connected with the issue of which rows are 
to be removed from T to obtain T, i.e. which nonzero elements of and M^> are left complex. 
We present some more theorems which answer this question. The general method is to construct 
all independent complex products of elements of M^' and which are rephasing-invariant with 
respect to (|5|), (|6|) and, if Xji?,' s are present, also ([?]). These must involve an even number of such 
elements, since for each index on an M, there must be a corresponding index on an M* in order to 
form a rephasing-invariant product. Since in general, by construction, these have arguments 7^ 0, ir, 
each one implies a constraint which is that the set of 2n elements which comprise it cannot be made 
simultaneously real by any rephasings or the charged lepton or neutrino fields. (In special cases, a 
symmetry of a model may render some of these real,.) In the neutrino sector, although the Dirac 
mass matrix conserves the global lepton number U(l), while the Majorana mass matrices violate 
this symmetry, it is possible and convenient to treat v c - R and XkR together, as was done in (||), the 
RHS of which can be written simply as (1/2) f^' + h.c. The rephasings in (||) and (0) 

can be represented in a similarly unified way as fjR — ► e^j fj R . We thus construct the rephasing 
invariants 

n 

P 2n;j 1 k 1 ,...j n k n ;a L = Yi M jak a M a L (j a )k a ( 21 ) 
a=l 

where / = £ and ai is an element of the permutation group S n , and 



P (0) = TT M (0) M (0) * C22) 

2n;ilj 2 vj2n-li2n;o- 11 ha-lha <TO'2o-l)o - (j2o) V ' 



0=1 

where a € S^n- (The difference in structure is due to (||).) 

At quartic order, 2n = 4, (|2^) and ((2^) both yield the same complex invariant, which (in the 
notation of (|2l|)) is 

p(f) = p(f) _ M (f) M (f) m (/) *M (/) * (23) 



for / = £, 0. Note the symmetries 



p(f) — pU) (HA) 

jiki,j2k2 32k2,jik\ \ I 
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pU) _ pit)* (OK) 

jikij2k2 jxki,jik\ V °) 



(whence also j 2 ^ 2 = Pj 2 l* j 1 k 2 ) f° r / = ^® and, for / = 0, the additional symmetry 

p(0) = p (o) f26) 

3132,3334 3231,3433 V ' 

In a shorthand notation, (23) for f = £ may be denoted £££*£*. It is useful to classify the complex 
p(0) invariants according to which submatrices — M^ L \ M^ R \ and/or — in they involve. 

We find five possible (complex) structures: LLL*L*, RRR*R*, LDL*D*, RDR*D*, and LRD*D*. 



Explicitly, the first two are given by (23) for / = L,R, respectively, and the others by the products 

n (L/) . =M {L) M U) M (L) *M {fnn) * (27) 

for f = D, 
and 

Q {LRDD) _ M (L) M (R) M (D)* -.(D)* / oq\ 

In terms of P^°^'s (and using (Q), these are 

(LD) _ (0) M (0) (0)* M (0)* _„(0) f3Q x 

ii iii2,j 3 mi - m h,j2 m j3,S+mi IVl ji,j 3 m j2,3+mi ~ r h 31, 33 3+mi 

- M (0) M (0) M (0) * M (0) * -P (0) f311 

q(lrdd) _ (o) M (o) M (o)* M (o)* _ „(o) f32l 

lL 3i32Mk 2 - 1V1 ji,j2 1V1 3+k 1 ,3+k2 1V1 h,3+ki 1V1 j2,3+k 2 - r j2 31, 3+fei 3+k 2 

Secondly, we construct a set of complex invariants connecting the Q = — 1 and Q = lepton 
sectors. At order 2n = 4 we find two independent types of complex invariants, which we denote as 

Q (Di) -m {D) M {£) M (D) *M (e) * (33) 

^3lkl,hmi ~ 1V1 jik 1 1VI j2m 1 M j 2 k 1 IVI jimi K 66 ) 

and Ilj^- 2 - 3mi as given by (|7j) for / = £. Note that 

(D£) = (D£)* , , 

^ilfclj2"ii ^32k\,jim 1 \ ) 

In the short notation, these are D£D*£* and L£L*£*. We also observe that 

u (Lf) = Q (L,f) , , 

At order 2n = 6, we find first the complex invariants f f f f* f* f* for / = £, 0, each of which 
involves only one of the two charge sectors. The various possibilities for the permutations in 
and (|2^) yield only one independent explicit form, which we denote simply by 

p U) _ M (f) M U) M U) m (/) *M (/) *M (/) * (36) 

r jiki,32k2,hk 3 ~ 1VI 3iki IV1 j2k2 1V1 j 3 k 3 IV1 j2k 1 IV1 j 3 k2 1VI j 1 k 3 
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The f = £ invariant is denoted ££££*£*£*. The / = invariants can be further characterized 
according to which submatrices of they involve; we find the following types: LLLL*L*L*, 
RRRR*R*R*, DDDD*D*D*, DDLD*D*L* , DDRD*D*R* , LLDL*L*D*, RRDR*R*D*, 
LRDL*R*D* , LLRL* D* D* , LRRD*D*R*, and LRDD*D*D*. (In our notation, ///'/*/*/'* 
and /'///'*/*/* are synonymous; also, LDDL*L*R* contains the same phase information as 
LLRL*D*D*, and so forth for the other conjugates of this type, which are hence not listed.) 
For any /, the invariants in (]36| ) satisfy 



P U) _ P U) __ p (f) ( o 7) 

jiki,j 2 k 2 ,j3k 3 jzkzjaksjiki jaksjikijzka \ ) 

and 

pU) _ p(/)* /oo\ 

Jlfcl ,32k2,33kz j3k2,j2kl,jlk3 V°°7 

For the cases where JVf"' is symmetric, i.e. M(V = M(°) or, for submatrices, , M( fl ) , we find 
the additional symmetry 

p(/) _ p(/) /on\ 

JX?2iJ3J4,J6J6 j2jl,j@j5,j4j3 \ ' 

We also find the following independent 6'th order complex invariants linking the charge Q = — 1 
and Q = sectors: ££L£*£*L*, LL£L*L*£\ ££D£*£*D*, DD£D*D*£*, £DL£*D*L*, £DR£*D*R\ 
£LR£* D* D* . The 6'th order mixed-charge invariants of the form ///'/*/*/*', and also the 6'th 
order Q = invariants which are of this form when written in terms of L, R and/or D submatrices, 
are given explicitly as 

Q(fff') _ M (f) M (f) M (f) M (f> M (f> M (f> (4Q) 

^hki,32k2,nm 1 ~ 1V1 hk 1 1V1 j 2 k 2 1V1 j 3 rn 1 1V1 j 1 k 2 1V1 j 3 k 1 M j 2 rn 1 l 4U J 

Some examples of explicit forms of other invariants include 

£DR£*D*R* : M^M^M^M^M^ (41) 
£LR£* D* D* : M^J^M^M^M^^ (42) 

£DL£*D*L* : M^M^M^M^M^M^ (43) 

M W M (D) M (L) M (i)* M (D)* (L)* 
M jiki m j2m 1 Ivl j 3 j4 m j3ki IVI j4m 1 lvI j 1 j2 ' erc ' 

Further discussion of higher order invariants is presented elsewhere []7j . 

Our theorems are as follows: For a given model, construct the maximal set of independent 
generally complex invariants; denote the number of these as Ni nv . Of the corresponding Ni nv 
arguments, a subset iVj a are linearly independent. Then (a) the number of unremovable, invariant, 
and hence physically meaningful phases N p is equal to Ni a ; (b) in the general case in which the Ni a 
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arguments are 7^ 0, tt, each implies a constraint that elements contained within the corresponding 
invariant cannot all be made simultaneously real; (b) this constitutes the complete set of constraints 
on which elements of and can be made simultaneously real. Often, the lowest-order 

nonvanishing complex invariants are quartic, but we give an exception below. 

In order to determine Ni a from the N^ nv independent complex invariants, we proceed as follows. 
Form the N eq vector of arguments of nonzero (and in general complex) elements of and , ^. 
For each complex invariant of a given order, X, arg(X) = J2fJ2j,k c ^/k ar 9(M( J? )jk, where the sum 
is over the N eq nonzero elements of and M^°\ These equations can be written as Z£ 



ILK 



where £ is the ^.^-dimensional vector £ = (arg(Xi), arg(Xj^ inv )) and Z is an Ni nv -row by 
A^ e9 -column matrix. Then rank(Z) = iVj . 

We illustrate our general theorems with some models^] We first consider the case n s = 0, i.e. 
no electroweak-singlet neutrinos. For models of this type, v is 6-dimensional, and T is N eq x 6 
dimensional. As an illustration, consider a model defined by 

/ E 12 

MW = E 2 i E 22 E 23 (44) 






E± 2 





E 2 i 


E 22 


E 2 3 





-E32 


E33 





L\ 2 







L 22 


L 2 3 





L 2 3 


L33 



M (i) = L 12 L 22 L 2 3 (45) 
y L 2 3 L33 J 

(We have verified that this model is able to fit current solar neutrino data and other established 
limits on neutrino masses and mixing.) For this model, there are N eq = 10 complex elements in 
and and corresponding rephasing equations. The matrix T is 10 x 6, and we find that 

it has rank 6. Our theorem ( |l9| ) then implies that there are N p = N eq — rank(T) = 4 unremovable 
phases in and M^ L \ We find Ni nv = 8 independent complex (quartic) invariants: P 22 33 = 

jp jp Tp* jp* p(0) p(^) j j r*2 tt(^) j jp j* jp* tt(^) r jp j* jp* 

^22-^33^32^231 ^22,33 _ ^22,33 _ -^22-^33-^23, ii 12,22 ~ -^12-^22 -^22-^12 , ii 12,32 ~ -^12-^32-^23^12 

n(^) j jp t* Tp* tt(^) j Tp t* ip* tt{L£) j 771 j* jp* tt(^) j jp j* jp* 

22,32 - L 22E 32 L 23 ^ 22 ll 22i 33 - -^22-C'33-^23- fc '23 ' ii 23,32 ~ ^23-^32-^33-^221 ii 23,33 ~ ^23-^33-^33^23 • 

From these we calculate the 8 x 10 dimensional Z matrix and find that it has rank 4, so that of 
the eight arguments of the complex invariants, there are Ni a = 4 linearly independent ones, in 
accord with the equality 2Vj a = N p and our result that N p = 4. Thus, for this model, it is not, in 
general, possible by any rephasings to make either the charged lepton or neutrino mass matrices 
real. Furthermore, the last complex invariant implies that it is not possible in general to make the 



Since our purpose is to illustrate our general theorems, we do not attempt to discuss an exhaustive set of viable 
models of neutrino masses and mixing. Indeed, given the unsettled status of atmospheric neutrino oscillations (claimed 
by some experiments, not by others), it is unclear what data a viable model must fit. Here we shall tentatively assume 
that the solar neutrino deficit does indicate neutrino oscillations but shall view it as premature to make any inference 
of new physics from the atmospheric neutrino oscillation data. 



8 



set of elements {L23 , L33 , E 23 , E33 } simultaneously real; the next-to-last implies the same for the 
set {L23, L33, E 22 , E32}, and so forth for each of the other complex invariants. 

A natural setting for models with electroweak-singlet neutrinos is an SO(10)grand unified theory 
(GUT), where a x'jL occurs, along with the SU(5) 10l and 5^ in the SO(10) 16l, for each of the 
three generations, so that n s = 3. An interesting model of this type was recently studied in a 
supersymmetric context in Ref. || (DHR). A pioneering study of a non-super symmetric SO(10) 
model which was similar but had more mass matrix parameters was carried out in Ref. || . In both 
cases the Higgs representations were restricted so as to obtain (complex) symmetric mass matrices. 
Here we shall consider two models. First, to exhibit the generality of our results, we will consider 
a generalization of the models of Refs. || and r in which M™' and are not assumed to be 

symmetric (recall from (|J) that is automatically symmetric) 

/ E l2 \ 

MW = E 2 i E 22 (46) 



M iR) 








E\2 





( 


E2\ 


E22 













£"33 


1 





D 12 







#21 





#23 


V 





D32 


D33 


/ 





Rl2 







R12 








V 








R33 



(47) 



(48) 



(Here all entries are in general complex.) Second, we shall consider the DHR model. We note 
that, to reduce the number of parameters, the Yukawa couplings were further restricted in the 
DHR model so that (besides the conditions AfW = M^ T and = M^) D 12 = -3R U 

and -D33 = — 3i?33. For the general model d46|)-(|4"8|) we calculate the 11 x 9 dimensional T matrix 
and find that it has rank 9, so that by our theorem (|l^ ) there are N p = N eq — rank(T) = 2 
unremovable, physically meaningful phases in the leptonic sector. Correspondingly, we find two 
independent complex invariants (which have independent arguments), so Ni nv = N{ a = 2. These 
are both of 6'th order: P 2 % 3fi6 = D 21 D 32 R 33 D 2 < 3 D* 33 R$2 and Q { 3 ^ 3 % = D 32 D 23 E 12 Dt 2 D* 33 E^ 2 . 
As can be checked explicitly here, making Afw and (complex) symmetric does not change the 

value of N p or make either of the complex invariants real. For the general model (and for the special 
case m which AfW and are symmetric), one thus has the constraints that the elements in the 

set {-D21, D 23 , -D32, -D33, R12, R33} cannot be made simultaneously real by any fermion rephasing, 
and similarly, the elements in the set {Z?i2, -D23, -D32, -D33, Ei 2 , E22} cannot be made simultaneously 
real by any rephasings; the first constraint shows that it is not possible to make both the Dirac 
neutrino mass matrix and the right-handed Majorana mass matrix simultaneously real 
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by any fermion field rephasings. For the light particle sector of the general model, our constraints 
allow one to place (i) both phases in the light neutrino mass matrix M v = M^M^~ l M^ T 
(with being made real); or (ii) one phase in M v and one in ; but (iii) not both in . In 
contrast, for the DHR model, it was noted || that there is only one (complex) phase in the leptonic 
sector (and indeed in the model as a whole). In terms of our formalism, this follows because, as 
a result of the relation between and given above, 53 66 is actually real, so that its 

phase, while still an unremovable, invariant, and hence meaningful phase, is no longer complex. It 
is still technically true that N p = Ni a = 2 in this model, but only one of these phases is complex. 
This model (in both its general form and the above specializations) also serves as an example of 
one in which there are no complex quartic invariants. We have shown that this is true of the quark 
sector of the model, where we found Q one 6'th order invariant, which, written in terms of Yukawa 
matrix elements, is $32*23 12 = ^3 ( 2 " ) ^2? )y i? )y i2 t) * y 33 t) * y 2 ( 2 i) *-0 Further results and applications will 
be reported elsewhere Q. 

This research was supported in part by NSF grant PHY-93-09888. 
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